Superconductor coupled to two Luttinger liquids 
as an entangler for electron spins 
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We consider an s-wave superconductor (SC) which is tunnel-coupled to two spatially separated 
Luttinger liquid (LL) leads. We demonstrate that such a setup acts as an entangler, i.e. it creates 
spin-singlets of two electrons which are spatially separated, thereby providing a source of electronic 
Einstein-Podolsky-Rosen pairs. We show that in the presence of a bias voltage, which is smaller 
than the energy gap in the SC, a stationary current of spin-entangled electrons can flow from the SC 
to the LL leads due to Andreev tunneling events. We discuss two competing transport channels for 
Cooper pairs to tunnel from the SC into the LL leads. On the one hand, the coherent tunneling of 
two electrons into the same LL lead is shown to be suppressed by strong LL correlations compared to 
single-electron tunneling into a LL. On the other hand, the tunneling of two spin-entangled electrons 
into different leads is suppressed by the initial spatial separation of the two electrons coming from 
the same Cooper pair. We show that the latter suppression depends crucially on the effective 
dimensionality of the SC. We identify a regime of experimental interest in which the separation of 
two spin-entangled electrons is favored. We determine the decay of the singlet state of two electrons 
injected into different leads caused by the LL correlations. Although the electron is not a proper 
quasiparticle of the LL, the spin information can still be transported via the spin density fluctuations 
produced by the injected spin-entangled electrons. 

PACS numbers: 73.40.Gk, 71.10.P, 74.10.P 
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I. INTRODUCTION 



Pairwise and non-local entangled quantum states, so- 
called Einstein-Podolsky-Rosen (EPR) pairs repre- 
sent the fundamental resource for quantum communica- 
tion schemes like dense coding, quantum teleportation 
or quantum key distribution [||, or more fundamentally, 
they can be used to test Bell's inequalities Experi- 
ments have tested Bell's inequalities [|j , dense coding , 
and quantum teleportation |j using photons, but to 
date no experiments for massive particles like electrons 
in a solid state environment exist. This is so because it is 
difficult to first produce entangled electrons and also to 
detect them afterwards in a controlled way due to other 
electrons interacting with the entangled pair. On the 
other hand, the spin of an electron was pointed out to 
be a most natural candidate for a quantum bit (qubit) 
H, 0. This idea was supported also by experiments 
which show unusually long dephasing times for electron 
spins in semiconductors (approaching microseconds) and 
phase coherent transport up to lOO^m ^ In ad- 
dition, the electron also possesses charge which makes it 
well suited for transporting the spin information ^ . 
Further, the Coulomb interaction between the electron 
charges can be exploited to spatially separate the spin- 
entangled electrons resulting in electronic EPR pairs. As 
first pointed out in Refs. |jl^, |l^, such electronic EPR 
pairs can be used for testing Bell inequalities and for 
quantum communication schemes in the solid state. The 
first step towards this goal is to have a scheme by which 
the electrons can be reliably entangled. One possibility is 
to use coupled quantum dots . Alternatively, we 



recently proposed an entangler device [ pl3| , which creates 
mobile and non-local spin-entangled electrons, consisting 
of an s-wave superconductor (SC), where the electrons 
are correlated in Cooper pairs with spin-singlet wavefunc- 
tions p^ . The SC is tunnel- coupled via two quantum 
dots in the Coulomb blockage regime to two spa- 
tially separated Fermi liquid leads. By applying a bias, 
a stationary current of spin-entangled electrons can flow 
from the SC to the leads. The quantum dots are used to 
mediate the necessary interaction between the two elec- 
trons initially forming a Cooper pair in the SC so that the 
two electrons tunnel preferably not into the same but in- 
stead into different leads. This entangler then satisfies all 
requirements to detect the entanglement via the current 
noise in a beam splitter setup [ |l8[ . It is straightforward 
to formulate spin measurements for testing Bell inequal- 
ities (it is most promising to measure spin via charge 
§, |l|, 1^). We refer to related work Q H ||, which 
makes also use of Andreev tunneling, but in a regime op- 
posite to the one considered in Ref. |l6[ and here, where 
the superconductor/normal interface is transparent and 
no Coulomb blockade nor strong correlations are present. 

In the present work we propose and discuss an al- 
ternative realization of an entangler which is based on 
strongly interacting one-dimensional wires which show 
Luttinger liquid (LL) behavior. In comparison to our 
earlier proposal with quantum dots |l^ , we replace now 
the Coulomb blockade behavior of the dots by strong 
correlations of the LL. Well-known examples for LL can- 
didates are carbon nanotubes . The low energy exci- 
tations of these LL are collective charge and spin modes 
rather than quasiparticles which resemble free electrons 
like they exist in a Fermi liquid. As a consequence. 
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the single-electron tunneling into a LL is suppressed by 
strong correlations. The question then arises quite natu- 
rally whether these strong correlations can even further 
suppress the coherent tunneling of two electrons into the 
same LL, as provided by a correlated two-particle tunnel- 
ing event (Andreev tunneling), so that the two electrons 
preferably separate and tunnel into different LL leads. It 
turns out that the answer is positive. To address this 
question we introduce a setup consisting of an s-wave 
SC which is weakly tunnel- coupled to the center (bulk) 
of two spatially separated one-dimensional wires 1,2 de- 
scribed as Luttingcr liquids, see Fig. 1,2. In this model 
we calculate the stationary current generated by the tun- 
neling of a singlet (spin-entangled electrons), transferred 
from the SC into two separate leads (non-local process) 
or into the same lead (local process), 1 or 2. We show 
that the ratio of these two competing current channels 
depends on the system parameters and that it can be 
made large in order to have the desired injection of the 
two electrons in two separate leads, where, again, the two 
spins, forming a singlet, are entangled in spin space while 
separated in orbital space and therefore represent an elec- 
tronic EPR pair. It is well-known that tunneling of single 




FIG. 1: A possible implementation of the entangler setup: 
Two quantum wires 1,2, described as infinitely long Luttinger 
liquids (LL), are deposited on top of an s-wave superconductor 
(SC) with chemical potential fis- The electrons of a Cooper 
pair can tunnel by means of an Andreev process from two 
points ri and r2 on the SC to the center (bulk) of the two 
quantum wires 1 and 2 resp. with tunneling amplitude to- 
The interaction between the leads is assumed to be negligible. 

electrons into LLs is suppressed compared to Fermi liq- 
uids due to strong many-body correlation. In addition, 
we find now that subsequent tunneling of a second elec- 
tron into the same LL is further suppressed, again in a 
characteristic interaction dependent power law, provided 
the applied voltage bias between the SC and the LL is 
much smaller than the energy gap A in the SC so that 
single-electron tunneling is suppressed. The two-particle 
tunneling event is strongly correlated within the uncer- 
tainty time h/A, characterizing the time-delay between 
subsequent tunneling events of the two electrons of the 
same Cooper pair. In other words, the second electron 
of a Cooper pair is incluenced by the existence of its pre- 
ceding partner electron already present in the LL. This 



effect can also be interpreted as a Coulomb blockade ef- 
fect, similar to what occurs in quantum dots attached 
to a SC [|l^ Similar Coulomb blockade effects oc- 
cur also in a mesoscopic chiral LL within a quantum dot 
coupled to macroscopic chiral LL edge-states in the frac- 
tional quantum Hall regime ||2^. There, the Coulomb 
blockade- like energy gap is quantized in units of the non- 
interacting energy level spacing of the quantum dot and 
its existence is therefore a finite size effect, whereas in 
the present case, we will see that the suppression comes 
from strong correlations in a two-particle tunneling event 
which is present even in an infinitely long LL as consid- 
ered here. On the other hand, if the two electrons of a 
Cooper pair tunnel to different leads, they will prefer- 
ably tunnel from different points from the SC ri and r2 , 
with distance 5r = ri — r2 due to the spatial separation 
of the leads, see Fig. 1,2. We find that the current is 
exponentially suppressed if the distance 6r exceeds the 
coherence length ^ of a Cooper pair on the SC. This lim- 
itation poses no severe experimental restriction since ^ is 
on the order of micrometers for usual s-wave materials, 
and 6r can be assumed to be on the order of nanometers. 
Still, a power law suppression cx l/{kpSr)'^ , with kp be- 
ing the Fermi wavevector in the SC, remains and is more 
relevant. We show, however, that in lower dimensions of 




FIG. 2; An alternative implementation of the proposed en- 
tangler set-up: The two quantum wires 1,2 with chemical po- 
tential fj,i , described as infinitely long Luttinger liquids (LL) , 
are tunnel-coupled with amplitude to from two points xi and 
X2 to two points ri and r2 of a superconducting (SC) tip with 
chemical potential fis- 



the SC this suppression is less pronounced (with smaller 
powers). Further, we then discuss the decay in time of 
a spin-singlet state injected into two LL, one electron in 
each lead, due to the interaction present in the LL and 
find a characteristic power law decay in time at zero tem- 
perature. Despite this decay of the singlet state, the spin 
information can still be transported through the LL wires 
via the spin-density fiuctuations created by the injected 
electrons. 
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II. MODEL AND HAMILTONIAN 

We consider an s-wave superconductor (SC) which is 
weakly tunnel-coupled to the center (bulk) of two spa- 
tially separated Luttinger liquid (LL) leads (see Fig. 1,2). 
The Hamiltonian of the whole system is represented as 
H = Ho + Ht with Ho = Hs + J2n=i.2 ^Ln describing 
the isolated SC and LL-leads 1,2 respectively. Tunneling 
between the SC and the leads is governed by the tun- 
neling Hamiltonian Ht- Each part of the system will be 
described in the following. 

The s-wave SC with chemical potential /^s is described 
by the BCS-Hamiltonian||| 

J?s-MsiV5-^i?k7L7ks, (1) 

where s = (t, i) and Ns = X^ks "^L'^ks is the number op- 
erator for electrons in the SC. The quasiparticle operators 
7ks describe excitations out of the BCS-groundstate |0)s 
defined by 7ks|0)5 = 0. They are related to the electron 
annihilation and creation operators Cks and cj^^ through 
the Bogoliubov transformation |2^ 

CkT = Wk7kT + "k7lki 
c-ki = '«k7-ki - Wk7kt ' 

where Wk = (1/V2)(1 -I- ^k/^^k)^/^ and Wk = (l/\/2)(l - 
^k/^'k)^^^ are the usual ECS coherence factors ||2^, and 
fk = Ck ~ Ms is the normal state single-elect ron energ y 
counted from the Fermi level /ig, and i?k = \/.^^ -I- is 
the quasiparticle energy. The field operator for an elec- 
tron with spin s is ^'^(r) — V^^^^^ X^k '^*'""'^ksj where V 
is the volume of the SC. 

The two leads 1,2 are supposed to be infinite one di- 
mensional interacting electron systems described by LL- 
theory. We only include forward scattering processes 
which describe scattering of electrons on the same branch 
(left or right movers). We neglect backscattering inter- 
actions which involve large momentum transfers of order 
2pp where is the Fermi wavevector in the LL [ p7[ . 
The LL-Hamiltonian for the low energy excitations of 
lead n = 1, 2 can then be written in a bosonized form as 

ii 

+L/2 

E / (^nL + ^(9.0„.)^), (3) 

where the fields n„(x) and (j)n{x) satisfy bosonic com- 
mutation relations [(l)ni^{x),Ilmii{x')] = i^nm^vn^ix — 
x'), and ^,1 is the chemical potential of the LL-leads 
(assumed to be identical for both leads), and Nn = 



j dxilil^g{x)tpns{x) is the number operator for elec- 
trons in LL n. The Hamiltonian (||) describes long- 
wavelength charge {v = p) and spin [v = a) density 
oscillations in the LL propagating with velocities Up and 
Ua, respectively. The velocities Ui, and the stiffness pa- 
rameters Kn depend on the interactions between the 
electrons in the LL. In the limit of vanishing backscat- 
tering, we have Mo- = vp and K„ = 1, and the LL 
is described by only two parameters Kp < 1 and Up. 
In a system with full translational invariance we have 
Up = vp/Kp. We decompose the field operator describ- 
ing electrons with spin s into a right and left moving 
part, ^nsix) = e'P^"Vn.+ (a;) + e-*P^"i^„,-(a;). The right 
(left) moving field operator ^ns+ (x) {ipns~ {x j) is then ex- 
pressed as an exponential of bosonic fields as E^, Is^ 



tpns±{x) = lim ^^i^expj ± ^(0„p(a;) 

+ S(t>na{x) T {Onp{x) "f sOna{x))^^, 

(4) 

where [4>niy{x),9rap.{x')] = -i(7r/2)5„,„5,.^sgn(x-a:') and 
therefore dxOnvix) = nlln^^x). The operators ?7±,„s are 
needed to ensure the correct fermionic anticommutation 
relations. In the thermodynamic (TD) limit (L — > oo), 
ri±,ns can be presented by Hermitian operators satisfying 
the anticommutation relation |p8| {77^ , rjr' } = 2Srr' , with 
r — ±,ns. We adopt the convention throughout the pa- 
per that s = -|-1 for s —'\, and s = —1 for s =|, if s has 
not the meaning of an operator index. 
Transfer of electrons from the SC to the LL-leads is de- 
scribed by the tunneling Hamiltonian Ht = X)n ^Tn + 
H.c, where Htu is defined as Htu = toJ2s "^ts"^ si'^n) ■ 
The field operator '^s{^n) annihilates an electron with 
spin s at point r„ on the SC, and i/'is creates it again 
with amplitude to at point x„ in the LL n which is nearest 
to r„, see Fig. 1,2. We assume that the spin is conserved 
during the tunneling process, and thus the tunneling am- 
plitudes to do not depend on spin, and, for simplicity, are 
the same for both leads n = 1,2. We remark that our 
point-contact approach for describing the electron trans- 
fer from the SC to the LL is the simplest possible descrip- 
tion but it captures presumably the relevant features of a 
real device. The scheme shown in Fig. 2 has a geometry 
which suggests that electrons tunnel from point r„ Xn 
which are closest to each other, due to the fact that to 
depends exponentially on the tunneling distance. In the 
setup shown in Fig. 1, a point-like tunnel contact between 
the SC and the LL might be induced by slightly bend- 
ing the quantum wires (e.g. nanotubes). If the contact 
area has a finite extension, we note that the two electrons 
preferably tunnel from the same point on the SC, when 
they tunnel into the same lead, since the two-particle 
tunneling event is coherent and shows a suppression in 
the probability already on a length scale given by 1/kp, 
as we discuss in detail below. 
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III. STATIONARY CURRENT FROM THE SC 
TO THE LL-LEADS 

We now calculate the current of singlets, i.e. pairwise 
spin-entangled electrons (Cooper pairs), from the SC to 
the LL-leads due to Andreev tunneling ||3l|, |3^ in first 
non-vanishing order, starting from a general T-matrix 
approach [ p3[ . We thereby distinguish two transport 
channels. First we calculate the current when two elec- 
trons tunnel from different points ri and r2 of the SC 
into different interacting LL-leads which are separated 
in space such that there is no inter-lead interaction. In 
this case the only correlation in the tunneling process is 
due to the superconducting pairing of electrons which 
results in a coherent two-electron tunneling process of 
opposite spins from different points ri and r2 of the 
SC, and with a delay time ~ h/A between the two 
tunneling events. Since the total spin is a conserved 
quantity [H, S^] = 0, the spin entanglement of a Cooper 
pair is transported to different LL-leads, thus leading to 
non-local spin-entanglement. On the other hand, if two 
electrons tunnel from the same point of the SC into the 
same LL-lead there is an additional correlation in the 
LL-lead itself due to the intra-lead interaction. It is the 
goal of this work to investigate how the transport current 
for tunneling of two electrons from the SC into the same 
LL-lead is effected by this additional correlation. 



IV. THE T-MATRIX 

We apply a T-matrix (transmission matrix) approach 
to calculate the current . The stationary current 
of two electrons passing from the SC to the leads is then 
given by 



I = 2eJ2Wf, 



(5) 



Here, Wfi is the transition rate from the superconductor 
to the leads, given by W/, = 27r|(/|r(e,)|i)|2j(e/ - Si) . 
Here, T{ei) — HTj-j^^zTni^i ~ ^o) is the on-shell trans- 
mission or T-matrix, with rj being a positive infinitesimal 
which we set to zero at the end of the calculation. The 
T-matrix can be expanded in a power series in the tun- 
neling Hamiltonian Ht, 



1 



Ei + iri- Hq 



HtT , (6) 



where Si is the energy of the initial state which, in 
our case, is the energy of a Cooper pair at the Fermi 
surface of the SC, = 2/1^. Finally, pi — {i\p\i) is 
the stationary occupation probability for the entire 
system to be in the state \i). We work in the regime 
A > /i > UbT, where p — ps ~ Pi the applied 
voltage bias between the SC and the leads, and T the 



temperature with ks the Boltzmann constant. The 
regime A > /i ensures that single electron tunneling 
from the SC to the leads is excluded and only tunneling 
of two coherent electrons of opposite spins is allowed. In 
the regime fi > ksT we only have transport from the 
SC to the leads, and not in the opposite direction. Since 
temperature is assumed to be the smallest energy scale 
in the system, we assume fc^T = in the calculation. 
The set of initial states \i), virtual states \v) and final 
states I/) consists of the BCS groundstate (GS) |0)5 
and excitations 7ks|0)s for the SC and a complete 
set |2^, |oj of energy eigenstates \Nnriy, {bnu}) of the 
LL-Hamiltonian H^n given in (H). iV„ri/ is the number 
of excess spin {v — a) and charge {v = p) in branch 
r relative to the state where all single-particle states 
are filled up to the chemical potential pi. The Bose 
operators bm, form a continuous spectrum describing 
collective spin and charge modes and will be introduced 
in and (|ll|). The GS of the LL is then |0, 0), which 
means that we have a filled Fermi sea and no bosonic 
excitations. The energy contribution of the excess 
charge and spin is included in the so-called zero mode 
(fc — 0) terms in the diagonalized Hamiltonian K^n ( p^ ) 
and are of no importance in the TD-limit (L — > oo) 
considered here, since the contribution of these terms 
due to an additional electron on top of the GS is 0{1/L) 
and is neglected in dl^. For a detailed description of 
the LL-Hamiltonian (pi) including the zero-modes see 
Appendix A. Since we want to calculate the transition 
rate for transport of a Cooper pair to the leads, the final 
states I/) of interest contain two additional electrons of 
opposite spins in the leads compared to the initial state 



V. CURRENT 7i FOR TUNNELING OF TWO 
ELECTRONS INTO DIFFERENT LEADS 

We first calculate the current for tunneling of two spin- 
entangled electrons into different leads. We expand the 
T-matrix to second order in Ht and go over to the inter- 
action representation by using S{e) '"^"^ 
and {v\{ei — Hq + iri)~^\v) = —i 

transforming the time dependent phases into a time de- 
pendence of the tunneling Hamiltonian we can integrate 
out all final and virtual states. The forward current Ii 
for tunneling of two electrons into different leads can then 
be written as 



{l/2n) j:^dte 



OO CXD CXD 



h = 2e lim^ J2 J "^^J"^^' J 6""^*'+*"^ 



-oo 



^i2t-t!-n, l^H^^Jt - t")Hl^,{t)HTn{t')HTn') , (7) 
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where (• • • ) denotes Trp{- - J^. The bias has been intro- 
duced in a standard way |34[| , and the time dependence 
of the operators in (R) is then governed by HTn{t) = 

giiKLr. + Ks)tjj^^^^-iiI^f„+Ks)t ^j^jj + Ks = Hlu + 

Hs — fJ-iNn — /is-^S- The transport process involves two 
electrons of different spins which suggests that the av- 
erage in is of the form (suppressing time variables), 

(•••) = J2ss'{H^Tm-s'HTm's'^Tn-sHTn's), whcrC Htus 

describes tunneling of spin s governed by Htu ■ The time 
sequence in (Q) contains the dynamics of the hopping of 
a Cooper pair from the SC to the LL-leads (one electron 
per lead) and back. The times t' and t" are delay times 
between subsequent hoppings of two electrons from the 
same Cooper pair, whereas t is the time between injecting 
and taking out a Cooper pair. We evaluate the thermal 
average in (|^) at zero temperature where the expectation 
value is to be taken in the groundstate of Kq which is 
the BCS groundstate of the SC and the bosonic vacuum 
of the LL-leads for a filled Fermi sea. We remark that 
since the interaction between the different subsystems 
{SC, Li, L2) is included in the tunneling-perturbation, 
the expectation value factorizes into a SC-part times a 
LL-part. In addition, the LL-correlation function factor- 
izes into two single-particle correlation functions due to 
the negligible interaction between the LL-leads 1, 2 (this 
will be not the case if two electrons tunnel into the same 
lead) . Note that in the TD-limit the time dynamics of all 
LL-correlation functions will be goverened by a Hamilto- 
nian that depends only on Bose operators (see (|l^)). The 
operators ?7±,rs commute with all Bose operators, and as 
a consequence r]±^rs are time independent. Therefore, in- 
teraction terms of the LL of the form tpa'4'i3'4'\'4's '^^^ 
written as 77077^77^775 x (Bose operators), where a,/?,7, <5 
are composite indices containing r = ±,ns. The correla- 
tion function in (Q) is then of the form 



= |<o|'' (V'«s(i-i")0(V'm-.(i-i')^i-.) 
X («'L,(r™, i - i")*I(r„, i)*-.(r™, i')*.(r»)) • (8) 



in (||) we then obtain 

1/2 (*t (r„, i _ t")*L,(r™, t)^'_,(r„„ t')*.(r„)) 



kk' 



(9) 



kk' 



where (5r = ri — r2 is the distance vector between the 
two tunneling points in the SC. The first sum in (^ de- 
scribes the (time-dependent) correlation of creating and 
annihilating a quasiparticle (with same spin), whereas 
the second term in (9|) describes correlation of creating 
two quasiparticles (with different spin). It is obvious 
that the second term describes processes which involve 
final states |/) in the T-matrix element {f\T{ei)\i) that 
contain two excitations in the SC and, therefore, does 
not describe an Andreev process. In the regime A > fj, 
such a process is not allowed by energy conservation. We 
will see this explicitly by calculating the integral over t 
which originates from the Fourier representation of the 
(5- function present in the rate Wfi . Similarly, for the cor- 
relator (^-L Jr„.,t-i")^'t(r„,<)^'_,(r,„,t')*.(r„)) in (|) 
we obtain (^ with a minus sign, and we have to replace 
t — t" hy t ~ t' — t", and t — t' hy t, in the second term 
of(|). 

To evaluate the LL-correlation functions in (g) we de- 
compose the phase fields (f>„^{x,t) and On^(x,t) into a 
sum over the spin and charge bosons (see also Appendix 
A), 



9„„{x,t) = sgn(p) 



2LK,\p\ 



^ipx ^-a\p\/2 



and 



-iUi, \p\t 



e 



iUi, \p\t\ 



(10) 



2L\p\ 

x(fe„.pe-'""l^l* + 6L-pe'""'^'*)- (11) 

The spin and charge bosons satisfy Bose-commutation 
relations, in particular [bni^p,bl^,^,p,] = Srr', where r = 
nup, and the LL-groundstate is defined as 5„i/p|0) ll = 0. 
The Hamiltonian (^) can then be written in terms of the 
6-operators as (see Appendix A) 



(12) 



up 



The 4-point correlation functions of the SC can be cal- 
culated by Fourier decomposing (r„, t) = V^^l"^ 

Ekl^ksTkse"*^"* -l-Wks7lk-se'^''*)e''"'": with Uks = "k, 
and = — Wkx = 'I'k- For the first correlation function discussed in Section IV and is explained in more detail in 



where we have subtracted the zero-point energy com- 
ing from the filled Dirac sea of negative-energy particle 
states. In all p-sums we will explicitly exclude 75 = as 
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Appendix A. To account for the p-dependence of the in- 
teraction, we apply a high momentum-transfer cut-off A 
on the order of 1/pf so that Ki,{p) = K^^ u^{p) = u^, 
for \p\ < 1/A and K„{p) — 1, Ui,{p) — vp for \p\ > 
1/A. By writing V«sr(x,t) = (2^a)-i/2^,,„, e'*— 
with r ~ ± and ^nsr defined according to (||), we 
can represent the single-particle LL-correlation func- 
tion as (V'nsr(a;,t)V'Lr> = (StTo) exp{ ($„sr (a;, t)^>„sr - 

5)72)} with the well-known result p3, 



nsr) 



1 K + i{vpt~rx) 

= — lim - 

27r a^o a -)- iivpt — rx) 

1 



n 



^ ^A + i{ui,t - rx) 



A2 



7^/2 

(13) 



where 7^ = {K^ + K~^)/A — 1/2 > is an interaction 
dependent parameter which describes the power law de- 
cay of the long time and long distance correlations. The 
factor in the first line of (|l^) is only important if one is in- 
terested in X, t satisfying \ vFt — rx\ < A and is a result of 
including the p-dependence of the interaction parameters 
K^j and it^ . The LL-correlation function has singular 
points as a function of time in the upper complex plane. 
It is now clear that the integration over t is only nonzero 
if the phase of e'"* in (0) is positive, i.e. w > 0. Since 
the phases of the terms containing (wk)^ in (^) depend 
also on i, the requirement for a nonzero contribution to 
the current from these terms requires 2/i — — E\^i > 0. 
This, however, is excluded in our regime of interest, since 
i?k + -Ek' > 2A, and therefore, we will not consider these 
terms any further in what follows. We remark that for 
r 7^ r' , the correlation function {ipnsr{x,t)il?^^^,) gives a 
negligible contribution in the TD-limit. This statement is 
true also for a finite size LL as long as the interaction pre- 
serves the total number of right- and left-movers. In our 
model the electrons tunnel into the same point Xn in LL 
n, i.e. a; = in (p^. In addition, LL 1 and 2 are assumed 
to be similar. We therefore have G^^^(x,t) = G^{t), and 
the current Ii can then be written as 



h = me\to\Vv') 



00 00 00 



lim / dt / dt' / dt"e-''(*'+*")e^(2t-t'-t")M 



X 

77^0 

-00 

X "k«kWk'«k' e-^^"*' e'E,,t"^rik+^')Sr 

kk' 

X {G\t - t")G\t - t') + G\t-t'- t")G\t)} . 

(14) 

We evaluate ( p^ ) in leading order in the small param- 
eter n/A and remark that the delay times t' and t" are 

restricted to t',t" ~ 1/A. This becomes clear if we set 



Sr = and express the contribution in ( p^ ) containing 
the dynamics of the SC as 



kk' 



(15) 



where H^^'' and iJQ^^* are Hankel functions of the first 
and second kind, and vs is the energy DOS per spin 
in the SC at /xg. For times t\t" > 1/A, the Hankel 
functions are rapidly oscillating, since for large x we 
have Ho (1/2) (2;) _ ^2/^ exp(±(ia; -_(7r/4))). In 
contrast, the time-dependent phase in (|l4| ) containing 
the bias fi suppresses the integrand in ( p^ ) only for times 
\2t -t' - t"\ > > 1/A. Being interested only in the 
leading order in /i/A, we can assume that \t\ > t',t" 
in the current formula (|lj), since the LL correlation 
functions are slowly decaying in time with the main 
contribution (in the integral) coming from large times 
\t\. In addition, since 1/fi > A/vp, we can neglect the 
term containing the Fermi velocity vp in (|l|). To test 
the validity of our approximations we first consider the 
non-interacting limit with K^, — 1 and u^, = vp, for 
which an analytic expression is also available for higher 
order terms in /i/A. 



VI. NON-INTERACTING LIMIT FOR 
CURRENT Ji 

Let us first consider a ID-Fermi liquid (i.e. K^, = 1 and 
Ui, = vp), and evaluate the integral over t in Eq. ( p^ in 
the non- interacting limit. The LL-correlation functions 
simplify to G^{t) = (l/27r) lima^o 1/(q: + ivpt), and we 
are left with the integral 



e^^tM {G\t- t")G\t - t') + G\t-t' - t")G\t)} 



1 



(a -I- ivF{t -t' - t")) {a + ivpt) 



(a + ivp{t - t")) {a + ivp{t - t')) 

(16) 



which can be evaluated by closing the integration contour 
in the upper complex plane. Inserting then the result into 
Eq. (|l|) we get 



00 00 



h = (32e|io|VV^^) lim^y dt' j di" e~''(*'+*") 


X "k^k^^k'i^k' e-'^"*' e^^^'*"e'(''+'^')^"^ 

kk' 



1 \sin{{t" -t')^i) sm{{t' +t")n) 



nvr 



t" - t' 



+ 



t' + t" 



(17) 
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The sine-functions in (|T^) can be expanded in powers of 
/Lt, and for /i < A it is sufficient to keep just the leading 
order term in fi since the integrals over t',t" have the 
form 



(77 ± i£;k)"+i 



(18) 



where n — 1,2,3 - • •. Since > A, higher powers in 
t',t" produce higher powers in /u/A, and, as expected, 
we can therefore ignore the dependence on t',t" in the 
LL-correlation functions. In contrast to this, when we 
consider the current for tunneling of two electrons into 
the same (interacting) LL-lead, we will see that the two- 
particle correlation function will not allow for such a sim- 
plification. In leading order in /i/A, the integrals over 
t',t" are evaluated according to (18) with n — 0, and we 
get an (effective) momentum-sum for the SC-correlations 
(^jj. ^^g^ cos (k • (5r))^. To evaluate this sum we use 
UkVk = A/(2i?k) and linearize the spectrum around the 
Fermi-level /is, since the Fermi energy in the SC ef > A, 
with Fermi wavevector kp. We then obtain {Sr denotes 
\Sr\) 



— - — cos [k ■ dr) = —i^s — ; — e 

k 



kpSr 



(19) 



In ( |l9|) we have introduced the coherence length of a 
Cooper pair in the SC, ^ = wf/ttA. We finally obtain 
/j*, the current /i in the non-interacting limit. 



sm{kpdr) 


' I 


' 2Sr\ 


kpSr 


exp \^ 





(20) 



Here we have defined 7 = 4TTi's'^i\tQ\'^/LV, which is 
the dimensionless conductance per spin to tunnel from 
the SC to the LL-leads. The non-interacting DOS of the 
LL per spin vi is given by vi — L/nvp. [We remark 
in passing that this result agrees with a T-matrix calcu- 
lation in the energy domain |3^. In this case we sum 
explicitly over the final states, given by a singlet |/) = 

(1/V2)[ajp|a2gx ~ '^ipi'^2qT]l*^' '^^^'^^ a-operators de- 
scribe electrons in a non-interacting ID Fermi-liquid. 
Note that triplet states are excluded as final states since 
our Hamiltonian H does not change the total spin.] We 
see that the current J° gets exponentially suppressed on 
the scale of 1^, if the tunneling of the two (coherent) elec- 
trons takes place from different points ri and r2 of the 
SC. For conventional s-wave SC the coherence length ^ is 
typically on the order of micrometers and therefore this 
poses not severe experimental restrictions. Thus, in the 
regime of interest 6r < ^, the suppression of the current 
/? is only polynomial, i.e. cx {l/kp6r)'^ . It was shown 
pO| that a superconductor on top of a two-dimensional 
electron gas (2DEG) can induce superconductivity (by 
the proximity effect) in the 2DEG with a finite order 



parameter. The 2DEG then becomes an effective two- 
dimensional (2D) SC. More recently, it was suggested 
that superconductivity should also be present in ropes 
of single- walled carbon nanotubes which are one- 
dimensional (ID) systems. It is therefore interesting to 
calculate (19) also in 2D and ID. In the case of a 2D 



SC we evaluate J2u '^^^^ cos (k • Sr) in leading order in 
dr/n^, and we find 



— - — COS (k • or) 

k(2D) 



^^s{MkpSr) + 2Y, 



J2uikpSr) 



(21) 



where Ju{x) denotes the Bessel function of order v. For 
large kpSr, we get J^{kpSr) ~ y^2/TrkpSr cos{kpSr — 
(1^71/2) — (7r/4)), which allows an approximation 
of the right-hand side of (|2l|) for large kpSr by 



{7r/2)iys^/2/TTkpSr cos{kpSr-{n/4)){l-{2/TT) In 2). This 
result is exact to leading order in an expansion in l/kpSr. 
So asymptotically, the current decays only cx 1/kpSr. 
For Sr — 0, the bracket on the right-hand side of (|2l| ) 
becomes 1 as in the 3D-case. 

In the case of a ID-SC we obtain 

^ cos (k • Sr) = %s cos{kpSr) e-^^'''^^\ (22) 
i/k 2 



E 

k(lD) 



where there are only oscillations and no decay of the 
Andreev amplitude (for Sr/ir^ < 1). We see that the 
suppression of the current due to a finite separation 
of the tunneling points on the SC can be reduced 
considerably (or even excluded completely) by going 
over to lower dimensional SCs. 



VII. CURRENT Ji INCLUDING 
INTERACTION. 



We now are ready to treat the interacting case. Having 
obtained confidence in our approximation schemes from 
the non-interacting case above, we can neglect now the 
t',t" dependence of the LL-correlation function appear- 
ing in (p^), valid in leading order in fi/ A.. In this limit 
the i-integral considerably simplifies to 



(27r)2 / ^^6^(2*-*'-*")/^ 

^ — 00 

X {G\t - t")GHt - t') + G\t-t' - t")G\t)} 



2A2(7p-I-7<t) 



,27. 



00 n.=p,.((A/w.)+«i) 



271.+1 ■ 



(23) 



An analytical expression for this integral is available , 
and given in Appendix B. The treatment of the remaining 
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integrals over t' ,t" and the calculation of the Andreev 
contribution is the same as in the non-interacting case 
and we obtain for the current /i, in leading order in /z/A 
and in the small parameters 2A/i/u,y, 



h 



VF 



r(27p + 2) Up 



27p 



(24) 



In (gj) we used K„ — 1 and = vp- The in- 
teraction suppresses the current considerably and the 
bias dependence has its characteristic non-linear form, 
Ii oc (/i)^''"''^^, with an interaction dependent exponent 
27p -|- 1. The parameter jp is the exponent for tunnel- 
ing into the bulk of a single LL, i.e. p{s) ~ le]^" [ p8[ , 
where p{e) is the single particle DOS. Note that the cur- 
rent Ii does not show a dependence on the correlation 
time 1/A, which is a measure for the time separation be- 
tween the two electron tunneling-events. This is so since 
the two partners of a Cooper pair tunnel to different LL 
leads with no interaction-induced correlations between 
the leads. 



VIII. CURRENT 1-2 FOR TUNNELING OF 
TWO ELECTRONS INTO THE SAME LEAD 



The main new feature for the case where two elec- 
trons, originating from an Andreev process, tunnel into 
the same lead, is now that the 4-point correlation func- 
tion of the LL no longer factorizes as was the case before 
when the two electrons tunnel into different leads (see 
(^). In addition the two electrons will tunnel into the 
lead prefarably from the same spatial point on the SC, 
i.e. Sr = 0. We denote by I2 the current for coherent 
transport of two electrons into the same lead, either lead 
1 or lead 2. It can be written in a similar way as Ii (see 
(^) with the difference that now we consider final states 
with two additional electrons (of opposite spin) in the 
same lead (either 1 or 2) compared to the initial state. 
We then obtain for I2, 



tunneling electrons. We obtain in this case 

E (^T„-s' (t - t")HLs' mTn-s{t')HTns) 

s,s' 

= \to\' E (^-(^ - i")^«-.(t)V'L.(i')V'L) 

s 

X {^livr^t - i")*^,(r„, i)*_,(r„, i')«'.(r„)) 

X (*L,(r„, t - t")*I(r„, m-s{rn, i')*.(r«)).(26) 

The 4-point correlation functions for the SC in ([2^ ) 
are the same as in (^) for the case when the two 
electrons tunnel into different leads, except that now 
dr = 0. The (normalized) 4-point correlation func- 
tions in (H) for the LL are Gl^,^{t,t' ,t") = (■0„rs(i - 

t")i^r.r'-sml,.sit'Wnrs)l{^nrs{t " i")^L.> 

{'4'nr'-sit — t')'ipnj.,_^), which Can be calculated using sim- 
ilar methods as described above for the single-particle 
correlation function. After some calculation we get 



n 



A - iui,t" 



K + iuu{t - t' - t") 



u—p.a 

' (A + iu^{t -t' ~ t")){K + iuat) 
(A + iup{t -t' - t")){A + iupt) 

[K-iupt"){K + iupt') 



r 








U- 





(27) 



where -y^^rr' = £.v{{l/K^) + rr' - (1 + rr'))/4 with 
^p/cr = ±1. The exponent ^vrr' is related to 71^, in- 
troduced in the single-particle correlation function ( p^ ) 
via ^vrr' = iulu for r = /, and ^^rr' = S.vi'^lu + 
1 — K^)/2 for r ^ r' . For the other sequence — 

{■^Pnr-sit - t")^nr's{t)^ir-s{t')^ir's)l{'^nr-s{t - t' - 
t")^ir-s){^nr's{t)^ir's)^ We obtaiu 



00 CXD 00 



/2 = 4elini^ j dt j dt' j dt" e~ 

-00 

i{2t-t'-t")tj. I rT-\ 



n(t'+t") 



(25) 

where we have used that the leads 1 and 2 are identi- 
cal which results in an additional factor of two. Again, 
the thermal average is to be taken at T = 0, and this 
groundstate expectation value factorizes into a SC-part 
times a LL-part. However, the LL-part does not factor- 
ize anymore due to strong correlations between the two 



n 



A - iu^t" 



A + iUjyt' 



u—p^a 

(A + iu„{t - t")){K + iu„{t - t')) 



A + iui,{t ~ t') 

l + rr' 



(A + iup{t - t"))(A + iup{t - t')) 
{K-iUpt"){K + iUpt') 



{K - iu„t"){K + iu„t') 



(28) 



We remark that contributions from other combinations 
of left- and right- movers, as indicated in (^) and 
(^), are negligible. A contribution like {ipnr'sit — 
t")ipnr~s{t)ipir'-si'^')'^irs)r^r' IS Only nou-zcro if Spin 
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exchange between right- and left-movers is possible, but 
this is a backscattering process which we explicitly ex- 



J 



elude. Using (25) together with 
expression for I2 (with Sr ~ 0) 



we obtain a formal 



00 00 00 



h = 4e 



^Tn^sAjtoT^ lim J dt J dt' J dt" e-^^''+'"^ e^^'''''-'"^'^ui'\t"A)il^^\t'A) 
-00 

X Yl {Gl{t,t',t")G\t-t")G\t-t')- Gl,{t,t',t")G\t-t'-t")G\t)). (29) 



b=±i 



In ( |29D the meaning of the summation index is 5 = +1 
for rr' = +1, and 6 = — 1 for rr' ~ —1. We proceed 
to evaluate the current I2 with = 1 and = vp, 
i.e. "/arr' = 0. Sincc 7prr' > 0, wc see from the first 
line in and (||) that for \t\ > A/up,t',t", the full 
4-point correlation function is suppressed by a factor 
~ [t't" /{t'^)Y'p^^' compared to its factorization approxi- 
mation. To calculate the current I2 we assume that the 
time scales A/up and A/vp are the smallest ones in the 
problem. The times A/up and A/vp are both on the order 
of the inverse Fermi energy in the LL, which is larger than 



the energy gap A and the bias /i. By applying the same 
arguments as in Section ^ for the current Ii , we approx- 
imate the current I2, assuming |i| > t',t" > A/vf,A/i 



which is accurate in leading order in the small parame- 
ters /i/A, AA/uj/ and A^l/u^. In this limit we obtain for 
the two-particle correlation functions Gl^,^ = — 

Up^'"'''' {t't"y>--' /{A+iuptf^"--' . The current h for tun- 
neling of two electrons into the same lead 1 or 2 then 
becomes (for 6r = 0) 



J 



^ / 27r^gA|^op \ _f^2_^^y f at' f dt"e-^^^''+'"Ut't'T' ii^o \t"ml^\t' A) 







1 f e*^'^* 



Again, the i nteg rals appearing in (|3C|) can be evaluated 
analytically p2[ with the results given in Appendix B. 
Note that accordingto the two-particle LL-correlation 
functions and (p8|), we find that the dynamics com- 
ing from the delay times t' and t" cannot be neglected 
anymore, as was done in We evaluate ( |30| ) in leading 
order in 2fiA/u,y and finally obtain for the current I2 



E 

fc=±i 



At, 



f) 



27pb 



(31) 



The interaction dependent constant Af, in (p^) is given 



by 

2^>^-i r(27p + 2) ,(lp, + l\ 

TT^ r(27pfc + 27p + 2)' I 2 j' ^^'^ 

which is decreasing for increasing the interactions in the 
LL leads and T{x) is the Gamma function. We remark 



that in ( pl| ) the current Ii is to be taken at (5r = 0. 
The non- interacting limit, I2 — Ii = Ii, is recovered by 
putting 7p = 7pf, = 0, and Up — vp. The result for I2 
shows that the unwanted injection of two electrons into 
the same lead is suppressed compared to Ji by a factor 
of A^{2fj,/ A)'^'^p+ if both electrons are injected into the 
same branch (left or right movers), or by A-{2fj,/ A)'^"'''- 
if the two electrons travel in different directions. Since 
it holds that 7p_ = 7p+ -|- (1 — Kp)/2 > jp+, it is more 
favorable that two electrons travel in the same direction 
than in opposite directions. The suppression of the cur- 
rent I2 by 1/A shows very nicely the two-particle corre- 
lation effect for the coherent tunneling of two electrons 
into the same lead. The larger A, the shorter is the delay 
time between the arrivals of the two partner electrons of 
a given Cooper pair, and, in turn, the more the second 
electron will be influenced by the presence of the first 
one already in the LL. By increasing the bias fj, the elec- 
trons can tunnel faster through the barrier due to more 
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channels becoming available into which the electron can 
tunnel, and therefore the effect of A is less pronounced. 
Also note that this correlation effect disappears when in- 
teractions are absent in the LL (7p = 7pf, = 0). 



IX. EFFICIENCY AND DISCUSSION 

We have established now that there exists indeed the 
suppression for tunneling of two spin-entangled electrons 
into the same LL-lead compared to the desired process 
where the two electrons tunnel into different leads. How- 
ever, we have to take into account that the process into 
different leads suffers also a suppression due to a finite 
tunneling separation Sr of the two electrons forming a 
Cooper pair in the SC. In Section Vl we showed that this 
suppression can be considerably reduced if one uses effec- 
tively low-dimensional SCs. To estimate the efficiency of 
the entangler we form the ratio /1//2 and demand that 
it is larger than one. This requirement is fulfilled if ap- 
proximately 



A, 



(33) 



where d is the dimension of the SC, and it is assumed 
that the coherence length ^ of the SC is large com- 
pared to Sr. The leading term of I2 is proportional to 
(2/i/A)^'''''+ describing the power-law suppression, with 
exponent 2jp+ = 27^, of the process where two elec- 
trons, entering the same lead, will propagate in the same 
direction. The exponent 7p4_ is the exponent for the 
single-particle tunneling-DOS from a metal (SC) into 
the center (bulk) of a LL. Experimentally accessible sys- 
tems which exhibit LL-behavior are metallic carbon nan- 
otubes. It was pointed out Q that the long range 
part of the Coulomb interaction, which is dominated by 
forward scattering events with small momentum trans- 
fer, can lead to LL behavior in carbon nanotubes with 
very small values of Kp 
imentally 



0.2 



0.3, as measured exper- 
45 1 and predicted theoretically |44 . This 
would correspond to an exponent 27p ^ 0.8 — 1.6, which 
seems very promising. In addition, single- wall nanotubes 
show similar tunneling exponents as derived here. The 
tunneling DOS for a single-wall nanotube is predicted to 
be /9(e) ~ ler' with 77 = {K-^+Kp-2)/8 Qg], which 
is half of 7p , and was measured H, H to be ~ 0.3 - 0.4. 
Similar values were found also in multiwall-nanotubes 
p6[ . It is known that the power-law suppression of the 
single-particle DOS is even larger if one considers tun- 
neling into the end of a LL. For single wall nanotubes 
one finds rjend ~ {Kp^ — l)/4 > 7/ Q, or for con- 
ventional LL-thcory again an enhancement by a factor 
of two 1^ . We therefore expect to get an even stronger 
suppression if the Cooper pairs tunnel into the end of 
the LLs. We remark that the non-locality of the two 
electrons could be probed via the Aharonov-Bohm oscil- 
lations in the current, when the leads 1,2 are formed into 
a loop enclosing a magnetic flux. Due to the different 



paths which the electrons can choose to go around the 
loop, we expect to see h/e and h/2e oscillation periods, 
as a function of magnetic flux, in the current like for non- 
interacting leads ||l6|] . Interference of contributions where 
the two electrons travel through different leads with con- 
tributions where they travel through the same lead then 
lead to the h/e oscillations, whereas interference of con- 
tributions where both electrons travel through the same 
arm 1 or 2 of the loop lead to the h/2e oscillations. The 
amplitudes of these oscillations must be related to the 
currents describing the interfering processes. We expect 
that the h/e oscillation contribution should be (x (I1I2)" 
and the h/2e oscillation contribution should be (x 
with an exponent a that has to be determined by ex- 
plicit calculations. In the non-interacting limit a should 
be 1/2 |l6|. The different periods then allow for an ex- 
perimental test of how successful the separation of the 
two electrons is. For instance, if the two electrons only 
can tunnel into the same lead, e.g. if kp5r is too large 
or the interaction in the leads too weak, then /i ~ and 
we would only see the h/2e oscillations in the current. 
The determination of the precise value of the exponents 
will be deferred to another publication since it requires 
a separate calculation including finite size properties of 
the LL along the hues discussed in 



X. DECAY OF THE ELECTRON-SINGLET DUE 
TO LL-INTERACTIONS 

We have shown in the preceding sections that the in- 
teraction in a LL-lead can help to separate two spin- 
entangled electrons so that the two electrons enter dif- 
ferent leads. A natural question then arises: what is the 
lifetime of a (non-local) spin-singlet state formed of two 
electrons which are injected into different LL-leads, one 
electron per lead? To address this issue we introduce the 
following correlation function 



F(r,t) = |(5(r,t)|5(0,0))r 



(34) 



This function is the probability density that a singlet 
state, injected at point r = {xi,X2) — and at time 
t = 0, is found at some later time t and at point r. 
Therefore, P(r, t) is a measure of how much of the initial 
singlet state remains after the two injected electrons have 
interacted with all the other electrons in the LL during 
the time interval t. Here 

\S{v,t)) = ^/^['ll)\^{xi,t)ll)l^{x2,t) 

- ^P\^{x^,t)^l^(x2,tm (35) 

is the electron singlet state created on top of the 
LL groundstates. The extra normalization factor 
V27ra is introduced to guarantee JdrP{r,t) = 1 
in the non-interacting limit and corresponds to the 
replacement of ipns by (/iTray^^ipns- The singlet- 
singlet correlation function factorizes into two single- 
particle Green's functions due to negligible interac- 
tion between the leads 1 and 2. Therefore we 
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have P(r,t) = (2™)2 n„ |(V^„.(a;„,t)V^L(0,0))| , with 
{^nsixn,t)^Pl{0,0)) = E.=± e''=-'^^"Gi,,(a:„,t). For 
simplicity we just study the slow spatial variations of 
\{^ns{xn,t)i'ns{0:^))\ ^nd obtain with (|l^ ) 



fluctuations can be expressed in a bosonic form (see Ap- 
pendix A) as 



Elim 
a-»0 



A'^ + {vpt - rxn)'^ 
o? + {vpt — rXn)^ 

A4 



n 



If we use tt6{x) = lima^o Q^/(a^ + 2;^) we can then 
write the remaining probability of the singlet as 

^('•'^)=n^E^W'^(^"-™^^*) (37) 

n r— ± 

with a time decaying weight factor of the ^-function 
Fit) 



n 



A2 



A2 + [VF - Utj)H'^ 

A" 




(A2 + {vFty^ + M^y + {2Auuty 



Without interaction we have F{t) — 1, which means 
that there is no decay of the singlet state. As interac- 
tions are turned on, we see that for times t > A/u^, the 
singlet state decays in time with approximately F{t) ~ 

n„_„„ ^ ^ f ) . This result to- 

gether with ( ^7|) shows that charge and spin of an electron 
propagate with velocity vp, whereas charge (spin) exci- 
tations of the LL propagate with Up {u„). In addition, 
we see that the probability P(r, t) shows an additional 
power-law decay ~ ( — / ^ y''" , with an interaction 

dependent exponent. We will show in the next section 
that although the singlet gets destroyed due to interac- 
tions, we still can observe charge and spin of the initial 
singlet via the spin and charge density fluctuations of the 
LL. 



XI. PROPAGATION OF CHARGE AND SPIN 

The charge and spin propagation as a function of time 
in a state can be described by the correlation func- 
tion {^\p{x,t)\^) for the charge, and {^\(7z{x,t)\'^) for 
the spin. The normal-ordered charge density operator 
for LL n is p„(a;„) = J2s '■ '^isixn)ipnsixn) := J2sr ■ 
i^nsrixn)i'nsr{xn) '■, if wc Only cousidcr the slow spa- 
tial variations of the density operator. Similarly, the 
normal-ordered spin density operator in z-direction is 
a^{xn) = J2sr^ ■ '4'lsrix„)tp„srixn) ■ ■ Thcse density 



PniXji) dx(l^np{Xji) 



and for the spin 



\ 7^/2 

(36) 



<^nix„) 



V2 



dx 



(39) 



(40) 



We now consider a state I'i') — ?/4isr(2;n)|0) where we 
inject an electron at time t — into branch r on top of 
the LL groundstate in lead n and calculate the time de- 
pendent charge and spin density fluctuations according to 
{0\ipnsr{xn)Pnix'n,t)'^isrixn)\0) for the charge and sim- 
ilar for the spin {0\^pnsr{xn)cr^{x'^,t)ipl^^^{xn)\0) . If we 
express the bosonic fieldoperators 4>nu and 9nu in terms 
of the boson modes shown in ([l^) and (|l^) and the Fermi 
operators according to the bosonization dictionary (jj) we 
obtain for the charge fluctuations 

nsr {Xn)Pn {x'n , t)''/'Lr iXn)\0) 

= i(l + rKp)S « - a;„ - Upt) 



+ -{1 - rKp)S (x'„ -Xn+ Upt) 



(41) 



and for the spin fluctuations 



nsr {•^n )^n i'^n ' )|0) 
= -(1 + rKa)5 {x'^ - Xn - Uat) 

+ -(1 - rK^)5 (x'„ - x„ + u^t) . 



(42) 



The results ( pdj ) and ( p^ ) are obtained by sending A — > 
and by using the same normalization convention for the 
electron operators as in (^5|). We see that in contrast to 
the singlet, the charge and spin density fluctuations in 
the LL created by the injected electron do not decay and 
show a pulse shape with no dispersion in time. This is 
due to the linear energy dispersion relation of the LL- 
model. In carbon nanotubes such a highly linear disper- 
sion relation is indeed realized, and, therefore, nanotubes 
should be well suited for spin transport. Another inter- 
esting effect that shows up in (^) and (^2|) is the different 
velocities of spin and charge, which is known as spin- 
charge separation. It would be interesting to test Bell 
inequalities Q via spin-spin correlation measurements 
between the two LL-leads and see if the initial entan- 
glement of the spin singlet is still observable in the spin 
density- fluctuations. Although detection of single spins 
with magnitudes on the order of electron spins has still 
not been achieved, magnetic resonance force microscopy 
(MFRM) seems to be very promising in doing so p8[ . 
Another scenario is to use the LL just as an intermediate 
medium which is needed to first separate the two elec- 
trons of a Cooper pair and then to take them (in general 
other electrons) out again into two (spatially separated) 
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Fermi liquid leads where the (possibly reduced) spin en- 
tanglement could be measured via the current noise in a 
beamsplitter experiment [|8|. Similarly, to test Bell in- 
equalities one can make then use of measuring spin via 
the charge of the electron |l^, ^ . In this context we 
finally note that the decay of the singlet state given by 
( ^ ) sets in almost immediately after the injection into 
the LLs (the time scale is approximately the inverse of 
the Fermi energy), but at least at zero temperature, the 
suppression is only polynomial in time, which suggests 
that some fraction of the singlet state can still be recov- 
ered. 



XII. CONCLUSIONS 



We proposed an s-wave superconductor (SC), coupled 
to two spatially separated Luttinger liquid (LL) leads, 
as an entangler for electron spins. We showed that the 
strong correlations present in the LL can be used to sep- 
arate two electrons, forming a spin-singlet state, which 
originate from an Andreev tunneling process of a Cooper 
pair from the SC to the leads. We have shown that the 
coherent tunneling of two electrons into the same lead 
is suppressed by a characteristic power law in the small 
parameter /i/A, where /i is the applied bias between the 
SC and the LL-leads, and A is the gap in the SC. On 
the other hand, when the two electrons tunnel into dif- 
ferent leads, the current is suppressed by the initial sep- 
aration of the two electrons. This suppression, however, 
can be considerably reduced by going over to effective 
lower-dimensional SC. We also addressed the question of 
how much of the initial singlet can be taken out of the LL 
at some later time, and we found that the probability is 
decreasing in time, again with a power-law (at zero tem- 
perature) . Nevertheless, the spin information can still be 
transported through the wires by means of the (proper) 
spin excitations of the LL. 

While preparing this manuscript we have learned of 
related and independent efforts by S. Vishveshwara et 
al.pgjl who consider a similar setup as proposed here 
thereby arriving at similar conclusions. 
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APPENDIX A: FINITE SIZE 
DIAGONALIZATION OF THE 
LL-HAMILTONIAN 

In this Appendix we derive the diagonalized form of 
the LL-Hamiltonian (|^) including terms of order 1/i, 
which describe integer charge and spin excitations. For 
simplicity, we consider only one LL and will therefore 
suppress the subscript n for the leads. We start with the 
exact bosonization dictionary for the Fermi-operator for 
electrons on branch r ~ ± Eg, pCfl , 



tprsix) = lim 



Ur 



: exp 



V27ra 

-I- ^ (<ppix) + s(f)aix) - r{9p{x) + s9a{x))^ |. 



72^ 



(Al) 



The ?7r.s-operator (often denoted as Klein factor) is uni- 
tary and decreases the number of electrons with spin 
s on branch r by one. This operator also ensures the 
correct anticommutation relations for ^j.g{x). The nor- 
mal ordered charge density operator is p{x) — X^sr ■ 
'4}'l,g{x)'4)rs{x) :, where : : measures the correspond- 
ing quantity relative to the groundstate with all sin- 
gle particle states filled up to the chemical potential 
/i;. The normal ordered spin density operator is de- 
fined by a^{x) — X^sr * ■ '4'ls{^)i'rs{x) : . In addi- 
tion, one can define (bare) current densitiy operators 
for charge jp = ^^j, ril}'l.g{x)il}rs{x), and for the spin 
jcr = X^sr '''^''i^rs{^)^rs{x), respectively. Note that the 
current density has not to be normal ordered since its 
groundstate expectation value vanishes. The normal or- 
dered product : 'tl)l^{x)'4)rs{x) : is calculated according 
to 

: ii)l,{x)iprs{x) := lim : i^^x + Ax)^rsix) : . (A2) 

Ax— )-0 

By expanding the operator product in ( [A2| ) within the 
normal-order sign, the right-hand side of ( |A2| ) equals 
{l/2TT)d,{(l)p{x) + scj),{x) - r{0p{x) + s9„{x)))/V2, from 
which one easily finds 

p{x) = —dx(l)p{x), a'^ix) = —da;(j)a{x), (A3) 

and for the current densities 

jp{x) = -V2np(a;), Mx) = -V2n,{x). (A4) 

The field Il,y{x) is related to 9^{x) by dxdu{x) — tt 11^(2;). 
We decompose the phase fields into 4'u{x) = (t>^{x) + 
(/)°(a;) and Il,y{x) — Uf^ {x) + 11° (x), where the part with 
non-zero momentum (j)^ and 11^ can be expanded in a 
series of normal modes 



p#0 



(A5) 
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and for the canonical momentum 




-"l^l/'(6.p-?>Lp)- (A6) 



These fields have to satisfy bosonic commutation rela- 
tions [(l)^{x),U^{x')] = iS,,^iS{x - x') - 1/L), which 

in turn demands [bup,bj^p,] = 5y^5pp' and [6^p,6^p'] = 
[6^p,6|^p,] = 0. The zero mode parts 4>1 and 11° can 
be found by considering the integrated charge (spin) 
and charge- (spin)-currents, respectively. For instance 
the integrated charge density X^rs — J dxpix) ~ 
(V2/7r)(0,(L/2) - 0p(-i/2)) = (x/2/7r)(0^(L/2) - 
(fPp{~L/2)). Similar results are obtained for the other 
density operators, which then implies the zero modes 
to be (t)l{x) = {t:/L){N+^ + 7V_,^)a; and Hi = 
-{1/L){N+, - iV_,,), where N^p/, = {N^^ ± Nri)/V2. 
The LL-Hamiltonian (||) is then diagonalized by the fol- 
lowing expansion of the bosonic fields 



(j)y{x) 



(A7) 



and for the canonical conjugate momentum operator, 



n.(x) 



bl 



-(iV+,-7V_,,), 



(A8) 



where we have used w- 
Kl = Hl — fJ,iN we then obtain 



pi^ — \p\l K^T^. For the operator 



p^O.u 

A,, 



(A9) 



In ( |Ag| ) we have subtracted the zero point energy 
(1/2) '^p-^Q ^ u,^\p\, which originates from an infinite filled 
Dirac sea of negative energy particle states in the LL- 
model. The zero modes in (A7) and (A8) give rise to 
contributions of order 1/L in the Hamiltonian ( [A9| ), and 
they are also responsible for a shift of the Fermi wavevec- 
tor Pi?, appearing in the Fermionic field operator ipnsix), 
by a contribution of the same order. Since we are only 
interested in the thermodynamic limit, we have neglected 
the zero mode contributions in explicit calculations in the 
main text. 



APPENDIX B: EXACT RESULTS FOR THE 
TIME INTEGRALS 



In this Ap pen dix we 
integrals in (|2^) and 
variable t appearing in 



rive the exact results for the time 
BOf). The integrals over the time 
and (30) have the form 



gi2Mt 27re"^^^ f2u)^+^'~^ 
'^^ ^ r(n\ L Q + R; 2A(u, 1 - u-')y) , (Bl) 

J \-U.a 



with Q^R > 1, but in general this integral is valid for 
Q, R satisfying Re{Q + R) > 1 (see [|| p. 345). The func- 
tion r(a;) in (Bl) is the Gammafunction and i_Fi(a; 7; z) 
is the confluent hypergeometric function given by 



a{a + 1) 



z 

iFi{a; 7; z) = 1 + -— -I- , -,^ 01 
7 1! 7(7 -I- 1) 2! 



+ 



rameters 2fiA/up and 2/xAAv- The integrals over the 
delay times t' and t" in (^0|) contain Hankelfunctions 
of the first and second kind which are linear combina- 
tions of Besselfunctions of the first and second kind, i.e. 
H^^/^^(iA) = Jo(tA) ± iYo{tA). The inegrals over t' and 



(B2) ^" ( pO| ) are therefore linear combinations of integrals 



of the form 



In the main text we considered only the leading order 
term of i-Fi since higher order terms are small by the pa- 
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lim / dt ( 



-r,t 



YQ{tA)f = lim 



-T{6 +1){A'' 



V 



2N|-i(5+l) 



Qs (v/V'rf 



(B3) 



and 



lim 



dt e 



Jo(tA)t*' - lim (r(5+l)(A2+,y2)-i(5+i)p^ (ry/v^fTA^)) 



(B4) 



This result is valid for S > -I (see [|2) p. 691). The 
functions Q and P are Legendre functions. The limit 
rj^O for Qsiv/VrFTA^) is (see || p. 959) 

QsiO) = -iv^sin(57r/2) (B5) 



and the limit 77 ^ for Ps{ti/ y/ff + A?) is 



^.(0) 



r(| + i)r(i^) 

5+1 
r(fti) 



^cos(<57r/2)i^V^- (B6) 
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